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The nuclear symmetry energy calculated in the RPA from the pairing plus sym- 
metry force Hamiltonian with equidistant single-nucleon levels is for mass number 
' yl = 48 approximately proportional to T(T-|- 1.025), where T is the isospin quantum 

. number. An isovector character of the pair field assumed to produce the observed 

\^ \ odd-even mass staggering is essential for this result. The RPA contribution to the 

CD ' symmetry energy cannot be simulated by adding to the Hartree-Fock-Bogolyubov 

energy a term proportional to the isospin variance in the Bogolyubov quasiparticle 
\ vacuum, and there are significant corrections to the approximation which consist 

' in adding half the isocranking angular velocity. The present calculation employs 

a smaller single-nucleon level spacing than used in a previous investigation of the 
model. 
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The so-called symmetry energy in the semi-empirical formula for nuclear 
masses is quadratic in the difference N — Z between the numbers of neutrons 
and protons [1]. In a more general view including isobaric analogue states, this 
may be seen as a quadratic dependence on the isospin quantum number T. 
Nuclei with small values of T have a larger binding energy than expected from 
a quadratic extrapolation from higher T- values [2]. In one parametrization of 
this extra binding energy, the so-called Wigner energy, a term linear in T is 
included in the total energy so that the symmetry energy becomes propor- 
tional to T(T + x) for some constant x. The best choice of x seems to be 
close to 1 [3,4], possibly with a preference for a slightly larger value [5]. Other 
parametrizations of the Wigner energy are employed in Refs. [2] and [6]. 
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The origin of the Wigner energy is a debated issue [7,8,2,9,10,11,12,13,14,15,6]. 
For a schematic, isobarically invariant Hamiltonian, I recently showed that the 
Hartree-Bogolyubov energy rises quadratically in T, whereas the T-dependence 
of the RPA correlation energy is dominated by a term equal to half the 
derivative of the Hartree-Bogolyubov energy with respect to T. The Hartree- 
Bogolyubov energy and the dominant T-dependent term in the RPA correla- 
tion energy thus give a symmetry energy proportional to T{T + 1) [16]. 

For the details of the theory, see Ref. [16]. The Hamiltonian is the sum of an 
independent-nucleon Hamiltonian with equidistant, fourfold degenerate single- 
nucleon levels, an isobarically invariant isovector pairing force and a symmetry 
force: 

fccrr 

Here, P is a pair annihilation isovector, T is the isospin, and G and k, are 
coupling constants. The single-nucleon energy takes Q equidistant values 
separated by r], and akar are single-nucleon annihilation operators, ka is k 
or k, where the bar denotes time reversal, and r distinguishes neutrons from 
protons. Hartree-Bogolyubov optimal states are derived from a Routhian 

R = H - XA^- fiT,, 

where Ay is the number of valence nucleons. The nucleon chemical potential 
A is placed midway between the lowest and the highest e^, so that (A^) ~ 2Q 
in the optimal state, and the isocranking angular velocity // varied to produce 
a range of {T^). 

The symmetry of even Ay/2 ± is imposed on the trial Bogolyubov quasi- 
particle vacuum. The model thus describes isobaric multiplets whose substate 
with Tz = T belongs to a doubly even nucleus. Due to this symmetry, the 
;z-coordinate of the gap isovector A = — G'(P) vanishes. By a suitable choice 
of phases, the optimal state is then an ordinary product of neutron and pro- 
ton BCS states with a common real neutron and proton gap A. For numeric 
convenience, A is kept constant with the variation of /x, which results in a 
neghgible variation of G. 

The RPA calculation employs the standard technique of a perturbative boson 
expansion of products of nucleon field operators [17]. By truncating the re- 
sulting expansions of H and R to second order, RPA operators H and R are 
obtained. The initial symmetries of H and R are recovered by replacing in H 
and R truncated expansions of Ay and T by formal variables which obey the 
mutual commutation relations of the exact operators and have the same com- 
mutators with other boson operators as the truncated ones. When {Ay) and 
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(Tz) are equal to eigenvalues of and respectively, the lowest eigenstate of 
the approximate Routhian thus obtained has = (A^) and T = = (Tz). 
This state is also an eigenstate of the approximate Hamiltonian with an eigen- 
value E that is the sum of two terms: the Hartree-Bogolyubov energy 

ho - {Ho} - + ^ 



and the RPA correlation energy 

i<j 



Here, AP — P — (P) and AT = T— (T), (Vi, denotes the eigenfrequencies of the 
boson system with Hamiltonian R, and b^j are the boson annihilation operators 
associated with pairs of Bogolyubov quasiparticles. Eq, E2 and E — Eq + E2 
are defined for any T = (T^) by these expressions. 

In the calculation of Ref. [16], which was intended to simulate the isobaric 
chain with mass number A — 48, I used the Fermi gas estimate of the single- 
nucleon level spacing f]. As pointed out by Satula and Wyss [18] and Glowacz, 
Satula and Wyss [5], this gives an unrealistically large value. I have there- 
fore repeated the calculation for two smaller values: r] = 4/(6A/(c MeV)/7r^) 
with c = 8 and 10 [19,20,21], which gives rj = 1.1 and 1.4 MeV respec- 
tively. Otherwise, the parameters of the present calculation are the same as 
in Ref. [16]: The number Q of equidistant, fourfold degenerate single-nucleon 
levels is equal to 24, the pair gap A = 12 MeV/\/A = 1.7 MeV [20], and 
the total coefficient of T^/2 in the Hartree-Bogolyubov part of the symmetry 
energy rj + K^ 2{134:AA-^ - 203.6^-^/^) MeV = 3.3 MeV [22]. 

The result for the smallest single-nucleon level spacing 1] = 1.1 MeV is shown 
in Fig. 1. A comparison with Fig. 1 of Ref. [16] shows that the change of rj 
does not alter the situation qualitatively. ^ Still, almost exactly £'0 — -E'o,t=o = 
{rj + k,)T'^/2. Since the pairing energy 2A^/G is essentially constant — it varies 
by less than . 1 MeV in the range of T considered — this means that even in the 
presence of pairing, the 'kinetic' contribution to the symmetry energy, that is, 
the contribution from {Hq), is in a very good approximation equal to riT^/2. 
This expression is easily shown to be exact for A = 0. 

Also like in the previous calculation, the RPA contribution to the symmetry 
energy E2 — £^2,t=o is dominated by the frequency of the normal mode 
associated with the direction of the isospin. This is indeed more accurately 



^ Due to a computation error, ///2 — (£'2 — -E2,T=o) as shown in Fig. 1 of Ref. [16] 
is too large by an almost constant factor about 1.2. 
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Fig. 1. £"0 — i?o,T=0) E2 — £'2,T=o and /x/2 as functions of T for r/ = 1.1 McV. 
Eq — Eq^t=o is scaled by a factor 1/3. 

true here than for larger values of ?]. Thus, from T = to 4 the sum of 
all the other terms in E2 only change by about .5 MeV for 77 = 1.1 MeV 
and about .8 MeV for rj = lA MeV as compared to about 1.5 MeV for the 
value f] = 2.1 MeV given by the Fermi gas estimate. The reason for this seems 
obvious: With the smaller level spacing 1], the isodeformation A/ rj is larger and 
the isorotation therefore more collective. This weakens the coupling between 
the collective and intrinsic degrees of freedom. It may be noticed that since 
A oc A^^/^ and 1] oc A~^, the isodeformation increases with A. Heavier nuclei 
should therefore show an even more collective isorotation. 

The identification of T with (Tz) provides an extension of E2{T) to negative T. 
Due to the isobaric invariancc of H , this results in an even function. The only 
term in E2 which is not analytic at T = is /i/2, which is replaced by — /x/2 
for negative T. It follows that E2 — -E2,t=o — A*/2 vanishes quadratically for 
T ^ 0+. Now, apart from a negligible correction due to the fact that A rather 
than G is kept constant in the calculation, we have ^ = (IEq/cIT = {rj + k)T. 
The curve in Fig. 1 is in fact well approximated by the second order polynomial 
E2 - ^2,T=o = {'n + «)^/2 - .033 MeV x T^. Altogether, we thus have in a 
good approximation 



E - Et=o = i^-^ .033 MeVj T{T + x) , 
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where 




.033 MeV 

iv + «)/2 



-1 



= 1.020. 



For 77 = 1.4 MeV, this analysis gives x — 1.031. These values are consistent 
with the possible indications in the data of x being slightly larger than 1 [5] . 

From my results in Ref. [16], Glowacz, Satula and Wyss infer x ^ .8 for 
T — 2 [5]. These authors may have calculated {E2 — E2^t=o)T/{Eq — Eq^t=o), 
which is evidently not a right way of extracting x from the energies. 

Nuclear masses are extensively compared with Hartree-Fock-Bogolyubov cal- 
culations. See for example Ref. [6] and references therein. If we neglect the 
change of optimal state by the exchange terms in the self-consistent single- 
nucleon potential and pair field, the difference between the Hartree-Fock- 
Bogolyubov and Hartrcc-Bogolyubov energies, the Fock term, is in the present 
theory precisely the first term in Expr. (1). This is a somewhat artificial term 
since it is cancelled by a part of the second sum in the second term. ^ In fact 
Eq. (1) may be recast into the form 



where Ei is the Bogolyubov quasiparticle energy. 

Satula and Wyss seem to suggest that the RPA correlation energy may be 
simulated by adding ?7(AT^)/2 to the Hartree-Fock-Bogolyubov energy [18]. 
Neglecting again the change of optimal state, this amounts in the present 
theory to replacing E2 by 



Since this is a modification of the Fock term, which is absent in Eq. (2), a 
relationship between E2 and E2 is not obviously anticipated. Fig. 2 shows 
that in fact they depend on T very differently. In particular, E'KT) is even 

and analytic at T = 0, so E2 — £^2 7 =0 vanishes quadratically for T ^ 0. Thus 
E2 does not give the empirical cusp at T = in the graph of the function 



^ In Ref. [16], in the third paragraph of the second colomn on page 289, 'second 
term in the second sum' should be as here 'second sum in the second term'. It is 
understood there that if A = 0, also G = 0. 




(2) 




E{T). 
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Fig. 2. {r]+G+K)T/2, E2-E2,t=o and E^-E* j.^(^ as functions of T for r/ = 1.1 MeV. 

The T-dependence of E2 mainly reflects that of the isospin variance (AT^). 
For A = 0, we have (AT^) = T. For A > 0, the diffuseness of the nucleon 
numbers makes a contribution to (AT^) which is largest for T = 0. (AT^) is 
an even function of T and for A > analytic at T = 0. This results in the 
hyperbola-like behaviour seen in the inset of Fig. 3 of Ref. [18] and reflected 
in Fig. 2. For A = 0, the first term in E2 has a T-dependent part GT/2, 
which arises from blocking of the neutron-proton pairing force when nucleons 
are promoted from proton to neutron levels. For A > 0, the asymptotic slope 
of the hyperbola gets a corresponding contribution from this term. Since G — 
.4 MeV <^ ri + K, — 3.3 MeV, this is small compared to the contribution of the 
second term. Evidently, this discussion of E2 applies analogously to the true 
Fock term without the term proportional to rj. 

As seen from Fig. 1, the RPA contribution to the symmetry energy is in the 
pairing plus symmetry force model largely taken into account by adding /i/2 to 
the self-consistent energy Eq. This corresponds to the approximation employed 
by Fraucndorf and Sheikh, who added fi/2 to an isocrankcd Hartree-Fock- 
Bogolyubov energy [23] . As pointed out by these authors, doing so is essentially 
equivalent to assigning the energy calculated for (T^) 



T+1/2 or ^r(T + l) 
to a state with isospin T, where T^, is some isospin coordinate. In fact, if 
Eo{{To,)) is the calculated energy, we have Eo [^TiT + 1) ) Eo{T + l/2) ^ 
Eo(T) + E'q(T)/2 = Eq(T) + /x/2. Such procedures are common in ordinary 
cranking calculations. Sec references in Ref. [23]. In a context of isocranking, 

[24] . In the pairing 



Satula and Wyss applied the constraint (T^) = yT{T + 1] 
plus symmetry force model, the T-dependence of the RPA correlation energy 
is seen to be taken into account by such recipes only in a crude approximation: 
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E2 has a T-dependence additional to that of ///2. It is in fact this additional 
T-dependence which makes x > 1. 

The presence of a strong isovector pair field is essential for getting x 1. 

For G = the Hartree-Fock approximation is in fact exact, and we have 
E - Et=o = ivT^ + kT{T + l))/2 = {r] + k)/2 X T{T + x), where x = 
K/{r] + K,) = .7 for 1] = 1.1 MeV. The leading contribution of the pairing 
force in perturbation theory {—GP^ ■ P) = G{T — 3fl)/2 adds a linear term 
GT/2 to the symmetry energy, but as long as G is below the threshold of pair 
condensation G ^/(2 + log (0/2)) it is much smaller than the term rjT/2 
required to get a; = 1. A pair gap A which is comparable to the single-nucleon 
level spacing rj is thus necessary in order to make the quasiparticle energies 
and the RPA frequencies except ii so independent of T that ii/2 becomes the 
dominant term in E2 — £'2,t=o- 

In conclusion, the pairing plus symmetry force model with equidistant single- 
nucleon levels introduced in Ref. [16] gives in the RPA for A = 48 a symmetry 
energy proportional to T{T + x), where x ~ 1.025. The large isovector pair 
field, whose strength is inferred from the observed odd-even mass staggering, 
is essential for this result. The RPA contribution to the symmetry energy is 
roughly equal to ^/2, where fi is the isocranking angular velocity, but there 
are significant corrections to this approximation. It is in fact these corrections 
which make a; > 1. The RPA contribution to the symmetry energy cannot be 
simulated by adding to the Hartree-Fock-Bogolyubov energy a term propor- 
tional to the isospin variance in the Bogolyubov quasiparticle vacuum. 
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